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ABSTRACT 

We consider the possibility that masses and gravitational potentials of galaxy cluster, 
estimated at X-ray wavelengths, could be explained without assuming huge amounts 
of dark matter, but in the context of /(i?)-gravity. Specifically, we take into account 
the weak field limit of such theories and show that the corrected gravitational potential 
allows to estimate the total mass of a sample of 12 clusters of galaxies. Results show 
that such a gravitational potential provides a fair fit to the mass of visible matter 
(i.e. gas + stars) estimated by X-ray observations, without the need of additional 
dark matter while the size of the clusters, as already observed at different scale for 
galaxies, strictly depends on the interaction lengths of the corrections to the Newtonian 
potential. 

Key words: alternative gravity, clusters of galaxies, dark matter. Clusters - X-rays 
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1 INTRODUCTION 

Since the pioneering work by (Zwicky 1933), the problems 
of high mass-to-light ratios of galaxy clusters and of the ro- 
tation curves of spiral galaxies have been faced by asking for 
huge amounts of unseen matter in the framework of Newto- 
nian theory of gravity. It is interesting to stress the fact that 
Zwicky addressed such an issue dealing with missing matter 
and not with dark matter. 

Later on, several versions of the the cold dark matter 
model (CDM) have been built starting from the assump- 
tion that a large amount of non-baryonic matter (i.e. matter 
non-interacting with the electromagnetic radiation) could 
account for the observations in the framework of the stan- 
dard Newtonian dynamics. Besides, by adding a further in- 
gredient, the cosmological constant A (Carroll et al. 1992; 
Sahni & Starobinski 2000), such a model (now ACDM) has 
become the new cosmological paradigm usually called the 
concordance model. 

In fact, high quality data coming from the measure- 
ments of cluster properties as the mass, the correlation func- 
tion and the evolution with redshift of their abundance (Eke 
et al. 1998; Viana et al. 2002; Bahcall et al. 2003; Bahcall 
& Bode 2003), the Hubble diagram of Type la Supernovae 
(Riess et al. 2004; Astier et al. 2006; Clocchiati et al. 2006), 
the optical surveys of large scale structure (Pope et al. 2005; 
Cole et al. 2005; Eisenstein et al. 2005), the anisotropics of 
the cosmic microwave background (de Bernardis et al. 2000; 
Spergel et al. 2003), the cosmic shear measured from weak 
lensing surveys (van Waerbecke et al. 2001; Refregier 2003) 
and the Lyman -a forest absorption (Croft et al. 1999; Mc- 



Donald et al. 2005) are evidences toward a spatially flat 
universe with a subcritical matter content and undergoing 
a phase of accelerated expansion. Interpreting all this infor- 
mation in a self-consistent model is the main task of modern 
cosmology and ACDM model provides a good fit to the most 
of the data (Tegmark et al. 2004; Seljak et al. 2005; Sanchez 
et al. 2006) giving a reliable snapshot of the today observed 
universe. 

Nevertheless, it is affected by serious theoretical short- 
comings that have motivated the search for alterna- 
tive candidates generically referred to as dark energy or 
quintessence. Such models range from scalar fields rolling 
down self interaction potentials to phantom fields, from phe- 
nomenological unified models of dark energy and dark mat- 
ter to alternative gravity theories (Peebles & Rathra 2003; 
Padmanabhan 2003; Copeland, Sami, and Tsujikawa 2006). 

Essentially, dark energy (or any alternative component) 
has to act as a negative pressure fiuid which gives rise to an 
overall acceleration of the Hubble fiuid. Despite of the clear 
mechanisms generating the observed cosmological dynamics, 
the nature and the fundamental properties of dark energy 
remain essentially unknown notwithstanding the great the- 
oretical efforts made up to now. 

The situation for dark matter is similar: its clustering 
and distribution properties are fairly well known at every 
scale but its nature is unknown, up to now, at a fundamental 
level. 

On the other hand, the need of unknown components as 
dark energy and dark matter could be considered nothing 
else but as a signal of the breakdown of Einstein Ceneral 



© 0000 RAS 



2 S. Capozziello, E. De Filippis, V. Salzano 



Relativity at astrophysical (galactic and extragalactic) and 
cosmological scales. 

In this context, Extended Theories of Gravity could be, 
in principle, an interesting alternative to explain cosmic ac- 
celeration without any dark energy and large scale struc- 
ture without any dark matter. In their simplest version, the 
Ricci curvature scalar R, linear in the Hilbert-Einstein ac- 
tion, could be replaced by a generic function f{R) whose true 
form could be "reconstructed" by the data. In fact, there is 
no a priori reason to consider the gravitational Lagrangian 
linear in the Ricci scalar while observations and experiments 
could contribute to define and constrain the "true" theory 
of gravity. For a discussion on this topic, sec (Capozziello 
& Francavigha 2008; Nojiri & Odintsov 2007; Sotiriou & 
Faraoni2008). 

In a cosmological framework, such an approach lead to 
modified Friedmann equations that can be formally writ- 
ten in the usual form by defining an effective curvature fluid 
giving rise to a negative pressure which drives the cosmic ac- 
celeration (Capozziello 2002; Capozziello, Carloni and Troisi 
2003; Carroll et al. 2004; Capozziello & Francaviglia 2008). 
Also referred to as f{R) theories, this approach is recently 
extensively studied both from the theoretical and the ob- 
servational point of view (see e.g. (Capozziello et al. 2003; 
Capozziello, Caxdone and Francaviglia 2006; Borowiec, God- 
lowski and Szydlowski 2006)). Moreover, this same approach 
has been also proposed, in early cosmology, as a mechanism 
giving rise to an inflationary era without the need of any 
inflaton field (Starobinsky 1980). 

All this amount of work has been, essentially, concen- 
trated on the cosmological applications and have convinc- 
ingly demonstrated that extended theories of gravity (in par- 
ticular f{R) gravity) are indeed able to explain the cosmic 
speed up and fairly fit the available data-sets and, hence, 
represents a viable alternative to the dark energy models 
(Capozziello, Cardonc and Troisi 2005). 

Changing the gravity sector could have consequences 
not only at cosmological scales, but also at the galactic and 
cluster scales so that it is mandatory to investigate the low 
energy limit of such theories. A strong debate is open with 
different results arguing in favor (Dick 2004; Sotiriou 2006; 
Cembranos 2006; Navarro & van Acoleyen 2005; AUemandi 
et al. 2005; Capozziello & Troisi 2005) or against (Dolgov & 
Kawasaki2003; Chiba 2003; Olmo 2005) such models. It is 
worth noting that, as a general result, higher order theories 
of gravity cause the gravitational potential to deviate from 
its Newtonian 1/r scaling (Stelle 1978; Kluske & Schmidt 
1996; Schmidt 2004; CUfton & Barrow 2005; Sobouti 2007; 
Mendoza & Rosas-Guevara 2007) even if such deviations 
may be vanishing. 

In (Capozziello, Cardone and Troisi 2007), the Newto- 
nian limit of power law f{R) = foR" theories has been in- 
vestigated, assuming that the metric in the low energy limit 
($/c'^ << 1) may be taken as Schwarzschild - like. It turns 
out that a power law term (r/rc)^ has to be added to the 
Newtonian 1/r term in order to get the correct gravitational 
potential. While the parameter f) may be expressed analyti- 
cally Eis a function of the slope n of the f{R) theory, Tc sets 
the scale where the correction term starts being significant. 
A particular range of values of n has been investigated so 
that the corrective term is an increasing function of the rar 
dius r thus causing an increase of the rotation curve with 



respect to the Newtonian one and offering the possibility to 
fit the galaxy rotation curves without the need of further 
dark matter components. 

A set of low surface brightness (LSB) galaxies with ex- 
tended and well measured rotation curves has been consid- 
ered (de Blok & Bosma 2002; de Blok 2005). These systems 
are supposed to be dark matter dominated, and successfully 
fitting data without dark matter is a strong evidence in fa^ 
vor of the approach (see also (Frigerio Martins & Salucci 
2007) for an independent analysis using another sample of 
galaxies). Combined with the hints coming from the cos- 
mological applications, one should have, in principle, the 
possibility to address both the dark energy and dark matter 
problems resorting to the same well motivated fundamen- 
tal theory (Capozziello, Cardone and Troisi 2006; Koivisto 
2007; Lobo 2008). Nevertheless, the simple power law }{R) 
gravity is nothing else but a toy-model which fail if one 
tries to achieve a comprehensive model for all the cosmolog- 
ical dynamics, ranging from the early universe, to the large 
scale structure up to the late accelerated era (Capozziello, 
Cardone and Troisi 2006; Koivisto 2007). 

A fundamental issue is related to clusters and superclus- 
tcrs of galaxies. Such structures, essentially, rule the large 
scale structure, and are the intermediate step between galax- 
ies and cosmology. As the galaxies, they appear dark-matter 
dominated but the distribution of dark matter component 
seems clustered and organized in a very different way with 
respect to galaxies. It seems that dark matter is ruled by 
the scale and also its fundamental nature could depend on 
the scale. For a comprehensive review see (Bahcall 1996). 

In the philosophy of Extended Theories of Gravity, the 
issue is now to reconstruct the mass profile of clusters with- 
out dark matter, i.e. to find out corrections to the Newton 
potential producing the same dynamics as dark matter but 
starting from a well motivated theory. This is the goal of 
this paper. 

As we will sec, the problem is very different with respect 
to that of galaxies and we need different corrections to the 
gravitational potential in order to consistently fit the cluster 
meiss profiles. 

The paper is organized as follows. In Section 2, also re- 
ferring to other results, we discuss the weak field limit of 
/(7?)-gravity showing that an f{R) theory is, in principle, 
quite fair to address the cluster problem (Lobo 2008). In 
Section 3, starting from the corrected potential previously 
derived for a point-like masses, we generalize the results to 
extended spherically symmetric systems as required for well 
shaped cluster models of the Bautz and Morgan classificar 
tion (Bautz 1970). Section 4 is devoted to the description of 
the general properties of galaxy clusters. In Section 5, the 
sample of galaxy clusters which we are going to fit is dis- 
cussed in details. Results are presented in Section 6, while 
Section 7 is devoted to the discussion and the conclusions. 



2 F(J^)-GRAVITY 

Let us consider the general action : 

A = j d''x^f^[}{R) + XCm\ , (1) 
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where f{R) is an analytic function of the Ricci scalar R, g 
is the determinant of the metric g'^j,, X = — - — is the cou- 
pling constant and Cm is the standard perfect-fluid matter 
Lagrangian. Such an action is the straightforward general- 
ization of the Hilbert-Einstein action of GR obtained for 
/(i?) = R. Since wc arc considering the metric approach, 
field equations are obtained by varying (1) with respect to 
the metric : 



1 X 



(2) 



where Tu,v = 



-2 



is the energy momentum tensor 



of matter, the prime indicates the derivative with respect to 
R and □ = ■a''^ . Wc adopt the signature (+, — , — , — ). 

As discussed in details in (Capozziello, Stabile and 
Troisi 2008a), we deal with the Newtonian and the post- 
Newtonian limit of f{R) - gravity on a spherically sym- 
metric background. Solutions for the field equations can be 
obtained by imposing the spherical symmetry (Capozziello, 
Stabile and Troisi 2007b): 



ds^ = goo{x'^ ,r)dx°^ + grr{x° ,'r)dr^ — r^dQ 



(3) 



where a;° = ct and dfl is the angular element. 

To develop the post-Newtonian limit of the theory, one 

can consider a perturbed metric with respect to a Minkowski 



background g^^y 
developed as: 



ViJ.i' +hfii,. The metric coefficients can be 



gu{t,r):^l+g^f>{t,r)+g[t\t,r) 
grr{t, r) ~ -1 -h gi-V {t, r) 
gee{t,r) = -r^ 



(4) 



where we put, for the sake of simplicity, c = 1 , a:: — ct —> t. 
We want to obtain the most general result without imposing 
particular forms for the /(i?)-Lagrangian. We only consider 
analytic Taylor expandable functions 



f{R) ~ /o + hR + /2ii' + faR^ + . 



(5) 



To obtain the post-Newtonian approximation of f{R) - grav- 
ity, one has to plug the expansions (4) and (5) into the field 
equations (2) and then expand the system up to the orders 
0(0), 0(2) and 0(4) . This approach provides general re- 
sults and specific (analytic) Lagrangians are selected by the 
coefficients fi in (5) (Capozziello, Stabile and Troisi 2008a). 

If we now consider the 0(2) - order of approximation, 
the field equations (2), in the vacuum case, results to be 



/irJi(2) - 2hg% + 8/2i?[? - firg'^ulr = 
2figiV - r[/irJi(2) 

-haSl - figS-l + Af2Rf + Af2rRfr] = 

I 2giV + r[2gi^l - rR^^^ + 2glV,r + rgS^,] = 

It is evident that the trace equation (the fourth in the sys- 
tem (6)), provides a differential equation with respect to the 
Ricci scalar which allows to solve the system at 0{2) - order. 
One obtains the general solution : 

(2) _ . 2GM 6i(t)e-^^ , ^2(t)e^^ 

9tt — 00 fT:;: ^ h -rr-prSTTT 



6(-C)3/2r 



2GM I '^l(t)l'V-t'-ij' 
fir 3(r 



2£t- 



where ^ 



6/2 



/i and /2 are the expansion coefficients 



obtained by the /(i?)-Taylor series. In the limit / — > i?, 
for a point-like source of mass M we recover the standard 
Schwarzschild solution. Let us notice that the integration 
constant So is dimensionless, while the two arbitrary time- 
functions Si{t) and 52 (i) have respectively the dimensions 
of lenght~^ and lenght~^\ ^ has the dimension lenght~'^ . 
As extensively discussed in (Capozziello, Stabile and Troisi 
2008a), the functions 5i{t) {i = 1,2) are completely arbitrary 
since the differential equation system (6) depends only on 
spatial derivatives. Besides, the integration constant So can 
be set to zero, as in the standard theory of potential, since 
it represents an unessential additive quantity. In order to 
obtain the physical prescription of the asymptotic ffatness 
at infinity, we can discard the Yukawa growing mode in (7) 
and then the metric is : 



ds 



2GM Siit)e 



dr 



1 + 



fir 3^r 

2GM _ Si{t){r^/^ + l)e-''^ 
fir 



3^r 



dr 



The Ricci scalar curvature is 



R 



Si{t)e- 



(8) 



(9) 



The solution can be given also in terms of gravitational 
potential. In particular, we have an explicit Newtonian-like 
term into the definition. The first of (7) provides the second 

order solution in term of the metric expansion (see the defini- 
tion (4)). In particular, it is gtt = 1 + 24^ 



= l + and 



then the gravitational potential of an analytic /(i?)-theory 
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GM _ (?i(t)e~''V^ 
fir 6^r 



(10) 



Among the possible analytic /(7?)-models, let us consider 
the Taylor expansion where the cosmological term (the 
above /o) and terms higher than second have been discarded. 
We rewrite the Lagrangian (5) as 



/(i?) ~ali^ + 02i^^ + .. 



(ll) 



and specify the above gravitational potential (10), generated 
by a point-like matter distribution, as: 



4>{r) 



where 



3GM / 
4air V 



(12) 



(13) 



L can be defined as the interaction length of the problem 
due to the correction to the Newtonian potential. We have 
changed the notation to remark that we are doing only a 

specific choice in the wide class of potentials (10), but the 
following considerations are completely general. 



3 EXTENDED SYSTEMS 

The gravitational potential (12) is a point-like one. Now we 
have to generalize this solution to extended systems. Let us 
describe galaxy clusters as spherically symmetric systems 
and then we have to extend the above considerations to this 
geometrical configuration. Wc simply consider the system 
composed by many infinitesimal mass elements dm each one 
contributing with a point-like gravitational potential. Then, 
summing up all terms, namely integrating them on a spher- 
ical volume, we obtain a suitable potential. Specifically, we 
have to solve the integral: 

$(r) = / r'^dr' / sin 9' dO' / dw' (t>{r') . (14) 
Jq Jo Jo 

The point- like potential (12)can be split in two terms. The 
Newtonian component is 

3GM 



<l>N{r) = 



(15) 



Aair 

The extended integral of such a part is well-known (apart 

3 ' 
4ai ' 

3 GM(<r) 



from the numerical constant expression. It is 



$jv(r) = 



4ai 



(16) 



where M(< r) is the mass enclosed in a sphere with radius 
r. The correction term 



(pdr) 



GM e-T 
4ai r 



(17) 



considering some analytical steps in the integration of the 
angular part, give the expression 



* Such a length is function of the series coefBcients, ai and 02, 
and it is not a free independent parameter in the following fit 
procedure. 



^c{r) = —-L 



f 

Jo 



dr'r' p{r') ■ 



— E — 



(18) 



The radial integral is numerically estimated once the mass 
density is given. We underline a fundamental difference be- 
tween such a term and the Newtonian one: while in the lat- 
ter, the matter outside the spherical shell of radius r does 
not contribute to the potential, in the former, external mat- 
ter takes part to the integration procedure. For this reason 
we split the corrective potential in two terms: 



if r' < r: 



^c,int{r)-- 



2nG 
4 

2nG 
4 



W-t'I \r+r'l 

■ L / dr'r' p(r') ■ 

Jo 

• L / dr'r' p(r') ■ e | — ^ — 1 

Jo \ J 



if r' > r: 



9c,ext{r)= ^ ■ L drrp{r)- 



27rG 



L dr r p{r ) ■ e 



-1 -I- 



4 / • ' ' \ r 

The total potential of the spherical mass distribution will be 

$(r) = $jv (r) + ^cint {r) + ^c,ext (r) (19) 

As we will show below, for our purpose, we need the gravi- 
tational potential derivative with respect to the variable r; 
the two derivatives may not be evaluated analytically so we 
estimate them numerically, once we have given an expres- 
sion for the total mass density p(r). While the Newtonian 
term gives the simple expression: 



d*N 
"dT 



(r) 



3 GM{<r) 
4ai 



(20) 



the internal and external derivatives of the corrective poten- 
tial terms are more involved. We do not give them explicitly 
for sake of brevity, but they are integral-functions of the 
form 



T{r, r') 



0(r) 



dr' f{r, 



(21) 



a(r) 



from which one has: 



dJ^jr, r') 
dr 



I3{r) 



dr'^i^ + 
dr 



-/(r,a(r))^(r) + /(r,/3(r))^(r) 



(22) 



Such an expression is numerically derived once the integra- 
tion extremes are given. A general consideration is in order 
at this point. Clearly, the Gauss theorem holds only for the 
Newtonian part since, for this term, the force law scales as 
1/r^. For the total potential (12), it does not hold anymore 
due to the correction. From a physical point of view, this is 
not a problem because the full conservation laws are deter- 
mined, for /(_R)-gravity, by the contracted Bianchi identities 
which assure the self-consistency. For a detailed discussion, 
see (Capozziello, Cardone and Troisi 2007; Capozziello & 
Francavigha 2008; Sotiriou & Faraoni2008). 
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4 CLUSTER MASS PROFILES 



Clusters of galaxies are generally considered self-bound grav- 
itational systems with spherical symmetry and in hydro- 
static equilibrium if virialized. The last two hypothesis are 
still widely used, despite of the fact that it has been widely 
proved that most clusters show more complex morphologies 
and/or signs of strong interactions or dynamical activity, es- 
pecially in their innermost regions (Chakrabarty et al. 2008; 
De Filippis et al. 2005). 

Under the hypothesis of spherical symmetry in hydrostatic 
equilibrium, the structure equation can be derived from the 
coUisionlcss Boltzmann equation 

d , / ^ 2n , 2pgas{r) 2 2 N / x d^ir) . . 

-^{pgas{r) 0-r-)+ \ {c^r-f^e.u,) = "Psa. W • (23) 

where 3> is the gravitational potential of the cluster, ar and 
(761, are the mass- weighted velocity dispersions in the radial 
and tangential directions, respectively, and p is gas mass- 
density. For an isotropic system, it is 

(Tr = cre,uj (24) 

The pressure profile can be related to these quantities by 

P(r) = a^pgasir) (25) 

Substituting Eqs. (24) and (25) into Eq. (23), we have, for 
an isotropic sphere. 



dP{r) 



-Pgas{r) 



d*(r) 



(26) 



dr "^""'^ ' dr 

For a gas sphere with temperature profile r(r), the velocity 
dispersion becomes 

kT{r) 



2 



(27) 



where k is the Boltzmann constant, p, w 0.609 is the mean 
mass particle and nip is the proton mass. Substituting Eqs. 
(25) and (27) into Eq. (26), we obtain 

d [ kT{r) \ d$ 

or, equivalently. 



d£ _ fcT(r) 

dr prripr 



d\npgas{r) ^ d\nT{r) 



(28) 



d In r ci In r 

Now the total gravitational potential of the cluster is: 
<E>(r) =$jv(r)-F$c(r) (29) 
with 

<E>c (r) = $c,mt (r) + <!>c,e.t (r) (30) 

It is worth underlining that if we consider only the stan- 
dard Newtonian potential, the total cluster mass Mci,n{j') 
is composed by gas mass -|- mass of gala^xies -|- cD-gal£ixy 
mass -|- dark matter and it is given by the expression: 

Mci,N{r)=Mgas{r) + Mgaiir) + Mcogaiir) + MoMir) 



kT( 



lirripG 



dill Pgasij-) fi hi T(r) 



d In r 



d hi r 



(31) 



Mci,N means the standard estimated Newtonian mass. Gen- 
erally the galaxy part contribution is considered negligible 
with respect to the other two components so we have: 



fcT(r) 
prup 



d\npgas{r) _^ d\nT{r) 



dlnr 



dlnr 



Since the gas-mass estimates are provided by X-ray obser- 
vations, the equilibrium equation can be used to derive the 
amount of dark matter present in a cluster of galaxies and 
its spatial distribution. 

Inserting the previously defined extended- corrected po- 
tential of Eq. (29) into Eq. (28), we obtain: 



d3>N 
""dT 



d^c 
dr 



kT(r) 
prUpr 



dlnpga.(r) ^ dlnT(r) 



dlnr 



dlnr 



(32) 



from which the extended- corrected mass estimate follows: 

4ai 2d$c . X 
M^,Ecir) + ^r — (r) = 



4ai 
'~3~ 



kT{r) 
lirUpG 



dlnpgasir) ^ d\nT{r) 



dlnr 



dlnr 



(33) 



Since the use of a corrected potential avoids, in principle, 
the additional requirement of dark matter, the total cluster 
mass, in this case, is given by: 

M^i,Ec{r) = Mgas{r) + Mgai{r) + McDgai{r) (34) 
and the mass density in the #c term is 

Pci,Ec{r) = Pgas{r) + pgal{r) + PCDgal{r) (35) 

with the density components derived from observations. 

In this work, we will use Eq. (33) to compare the bary- 
onic mass profile Md^Ecir), estimated from observations, 
with the theoretical deviation from the Newtonian gravi- 

4ai 2 d$c 

tational potential, given by the expression r (r). 

3G dr 

Our goal is to reproduce the observed mass profiles for a 
sample of galaxy clusters. 



5 GALAXY CLUSTER SAMPLE 

The formalism described in § 4 can be applied to a sample of 
12 galaxy clusters. We shall use the cluster sample studied 
in (Vikhlinin et al.2005; Vikhlinin et al.2006) which consists 
of 13 low-redshift clusters spanning a temperature range 
0.7 9.0 keV derived from high quality Chandra archival 
data. In all these clusters, the surface brightness and the 
gas temperature profiles are measured out to large radii, so 
that mass estimates can be extended up to rsoo or beyond. 



5.1 Gas Density Model 

The gas density distribution of the clusters in the sample 
is described by the analytic model proposed in (Vikhlinin 
et al.2006). Such a model modifies the classical /3— model 
to represent the characteristic properties of the observed X- 
ray surface brightness profiles, i.e. the power-law-type cusps 
of gas density in the cluster center, instead of a fiat core 
and the steepening of the brightness profiles at large radii. 
Eventually, a second /3— model, with a small core radius, is 
added to improve the model close to the cluster cores. The 
analytical form for the particle emission is given by: 



2 

■ no 



('■/'v) 



i 



(l-|-rVr§)3/3-«/2 {l + r-'/riy/'T 
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2 
"02 



(1 + rVr22)3/32 



(36) 



which can be easily converted to a mass density using the 

relation: 



= riT ■ firUp 



1.4 
12' 



(37) 



where nr is the total number density of particles in the gas. 
The resulting model has a large number of parameters, some 
of which do not have a direct physical interpretation. While 
this can often be inappropriate and computationally incon- 
venient, it suits well our case, where the main requirement 
is a detailed qualitative description of the cluster profiles. 
In (Vikhlinin et al.2006), Eq. (36) is applied to a restricted 
range of distances from the cluster center, i.e. between an 
inner cutoff rmin, chosen to exclude the central tempera- 
ture bin (« 10 ^ 20 kpc) where the ICM is hkely to be 
multi-phase, and rdet, where the X-ray surface brightness 
is at least 3cr significant. We have extrapolated the above 
function to values outside this restricted range using the 
following criteria: 

• for r < r-min, we have performed a linear extrapolation 
of the first three terms out to r = kpc; 

• for r > rdet , we have performed a linear extrapolation of 
the last three terms out to a distance f for which pgas{r) = 
Pc Pc being the critical density of the Universe at the cluster 
redshift: Pc = Pc,o ■ (1 + z)^. For radii larger than f, the gas 
density is assumed constant at pgas{r)- 

We point out that, in Table 1, the radius limit rvn,,, is almost 
the same as given in the previous definition. When the value 
given by (Vikhlinin et al.2006) is less than the cD-galaxy ra- 
dius, which is defined in the next section, we choose this last 
one as the lower limit. On the contrary, Vmax is quite differ- 
ent fi:om rdet'- it is fixed by considering the higher value of 
temperature profile and not by imaging methods. 
We then compute the gas mass Mgasij) and the total mass 
Mci,N{r), respectively, for all clusters in our sample, sub- 
stituting Eq. (36) into Eqs. (37) and (31), respectively; the 
gas temperature profile is described in details in § 5.2. The 
resulting mass values, estimated at r = rmax, are listed in 
Table 1. 



5.2 Temperature Profiles 

As stressed in § 5.1, for the purpose of this work, we need 
an accurate qualitative description of the radial behavior of 
the gas properties. Standard isothermal or polytropic mod- 
els, or even the more complex one proposed in (Vikhlinin et 
al.2006), do not provide a good description of the data at all 
radii and for all clusters in the present sample. We hence de- 
scribe the gas temperature profiles using the straightforward 
X-ray spectral analysis results, without the introduction of 
any analytic model. 

X-ray spectral values have been provided by A. Vikhlinin 
(private communication). A detailed description of the rela- 
tive spectral analysis can be found in (Vikhlinin et al.2005). 



5.3 Galaxy Distribution Model 

The galaxy density can be modelled as proposed by (Balicall 
1996). Even if the galaxy distribution is a point-distribution 



1 . X 10^ 



1 . X 10^ 



1 . X 10^ 



100 150 200 300 500 700 1000 1500 
R ( Kpc ) 

Figure 1. Matter components for A478: total Newtonian dynam- 
ical mass (continue line); gas mass (dashed line); galactic mass 
(dotted-dashed line); cD-galaxy mass (dotted line). 



instead of a continuous function, assuming that galaxies are 
in equilibrium with gas, we can use a /3— model, oc r~^, for 
r < Rc from the cluster center, and a steeper one, oc r~^'^, 
for r > Rc, where Rc is the cluster core radius (its value is 
taken from Vikhlinin 2006). Its final expression is: 



Pgal{r) 



Pgal,l • 

< 




Pgal,2 • 

V 





r<Rc 
r> Rc 



(38) 



where the constants Pgai,i and Pgai,2 are chosen in the fol- 
lowing way: 

• (Bahcall 1996) provides the central number density of 
galaxies in rich compact clusters for galaxies located within a 

1.5 h^^Mpc radius from the cluster center and brighter than 
ms -|- 2™ (where ma is the magnitude of the third brightest 
galaxy): ngai,o ~ lO'Vi"* galaxies Mpc^^. Then we fix Pgai,i 
in the range ~ 10"''' 10"^'' kg/kpc**. For any cluster obeying 
the condition chosen for the mass ratio gal-to-gas, we assume 
a typical elliptical and cD galaxy mass in the range 10^^ 
IO^^Mq. 

• the constant Pgal,2 has been fixed with the only require- 
ment that the galaxy density function has to be continuous 

at Rc- 

We have tested the effect of varying galax;y density in the 
above range ~ lO^** 10^'' kg/kpc^ on the cluster with the 
lowest mass, namely A262. In this case, we would expect 
great variations with respect to other clusters; the result is 
that the contribution due to galaxies and cD-galaxy gives a 
variation < 1% to the final estimate of fit parameters. 
The cD galaxy density has been modelled as described in 
(Schmidt 2006); they use a Jaffe model of the form: 



PCDgal = 



PO,J 



^ 1 + ^ 



(39) 



where rc is the core radius while the central density is ob- 
4 3 

tained from Mj = —'kRcPo,j- The mass of the cD galaxy 

has been fixed at 1.14 x 10^^ Afg, with Tc = -Re/0.76, with 
Re = 25 kpc being the effective radius of the galaxy. The 
central galaxy for each cluster in the sample is assumed to 
have approximately this stellar mass. 
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Table 1. Column 1: Cluster name. Column2: Richness. Column 3: cluster total mass. Column 4: gas mass. 
Column 5: galaxy mass. Column 6: cD-galaxy mass. Gas and total mass values are estimated at r = rmax- 

Coluiiiii 7: ratio of total galaxA' mass to ,y;as mass. Colum^i 8: minimum radius. Column 9: maximum radius. 



name 


R 


Ma, 


AT 


Mgaa 




M^Ug 


al 


gal 
gas 




f'max 






(Mq) 


( !l/f 
(™0 


) 


(Me) 


(Mq) 


(kpc) 


(kpc) 


A133 





4.35874 


•lO" 


2.73866 


10" 


5.20269 ■ 


1012 


1.10568 ■ 


1012 


0.23 


86 


1060 


A262 





4.45081 


•10" 


2.76659 


1012 


1.71305 ■ 


10" 


5.16382 ■ 


1012 


0.25 


61 


316 


A383 


2 


2.79785 


■10" 


2.82467 


10" 


5.88048 ■ 


1012 


1.09217 ■ 


1012 


0.25 


52 


751 


A478 


2 


8.51832 


•10" 


1.05583 


10" 


2.15567 ■ 


10" 


1.67513 ■ 


1012 


0.22 


59 


1580 


A907 


1 


4.87657 


10" 


6.38070 


10" 


1.34129 ■ 


IQis 


1.66533 • 


10" 


0.24 


563 


1226 


A1413 


3 


1.09598 


10" 


9.32466 


10" 


2.30728 ■ 


IQis 


1.67345 • 


10" 


0.26 


57 


1506 


A179r5 


2 


5.44761 


■10" 


5.56245 


10" 


4.23211 ■ 


10" 


1.93957 ■ 


10" 


0.11 


79 


1151 


A1991 


1 


1.24313 


10" 


1.00530 


10" 


1.24608 ■ 


10" 


1.08241 • 


10" 


0.23 


55 


618 


A2029 


2 


8.92392 


10" 


1.24129 


10" 


3.21543 ■ 


lOi-^ 


1.11921 ■ 


10" 


0.27 


62 


1771 


A2390 


1 


2.09710 


■10" 


2.15726 


10" 


4.91580 ■ 


10" 


1.12141 ■ 


10" 


0.23 


83 


1984 


MKW4 




4.69503 


■10" 


2.83207 


1012 


1.71153 ■ 


10" 


5.29855 ■ 


10" 


0.25 


60 


434 


RXJ1159 




8.97997 


■10" 


4.33256 


1012 


7.34414 ■ 


10" 


5.38799 ■ 


10" 


0.29 


64 


568 



We have assumed that the total gala^xy-component mass 
(galaxies plus cD galaxy masses) is Ri 20 25% of the gas 
mass: in (Schindlor 2004), the moan fraction of gas versus 
the total mass (with dark matter) for a cluster is estimated 
to be 15^20%, while the same quantity for galaxies is 3-7-5%. 
This means that the relative mean mass ratio gal-to-gas in a 
cluster is w 20-^25%. We have varied the parameters Pgai,i, 
pgai,2 and Mj in their previous defined ranges to obtain a 
mass ratio between total galaxy mass and total gas mass 
which lies in this range. Resulting galaxy mass values and 
sal 

ratios , estimated at r = rmax, are listed in Table 1. 

gas 

In Fig. (1), we show how each component is spatially 
distributed. The CD-galaxy is dominant with respect to the 
other galaxies only in the inner region (below 100 kpc). As 
already stated in § 5.1, cluster innermost regions have been 
excluded from our analysis and so the contribution due to 
the cD-galaxy is practically negligible in our analysis. The 
gas is, as a consequence, clearly the dominant visible com- 
ponent, starting from innermost regions out to large radii, 
being galaxy mass only 20 25% of gas mass. A similar be- 
havior is shown by all the clusters considered in our sample. 



5.4 Uncertainties on mass profiles 

Uncertainties on the cluster total mass profiles have been 
estimated performing Monte-Carlo simulations (Neumann 
1995). We proceed to simulate temperature profiles and 
choose random radius-temperature values couples for each 
bin which we have in our temperature data given by 
(Viklilinin ct al.2005). Random temperature values have 
been extracted from a Gaussian distribution centered on 
the spectral values, and with a dispersion fixed to its 68% 
confidence level. For the radius, we choose a random value 
inside each bin. We have performed 2000 simulations for 
each cluster and perform two cuts on the simulated profile. 
First, we exclude those profiles that give an unphysical neg- 
ative estimate of the mass: this is possible when our simu- 
lated couples of quantities give rise to too high temperature- 
gradient. After this cut, we have « 1500 simulations for any 
cluster. Then we have ordered the resulting mass values for 
increasing radius values. Extreme mass estimates (outside 



the 10 90% range) are excluded from the obtained distri- 
bution, in order to avoid other high mass gradients which 

give rise to masses too different from real data. The resulting 
limits provide the errors on the total mass. Uncertainties on 
the elcctrou-dcnKity profiles have not been included in the 
simulations, being them negligible with respect to those of 
the gas-temperature profiles. 



5.5 Fitting the mass profiles 

In the above sections, we have shown that, with the aid of 
X-ray observations, modelling theoretically the galaxy dis- 
tribution and using Eq. (33), we obtain an estimate of the 
baryonic content of clusters. 

We have hence performed a best-fit analysis of the theoret- 
ical Eq. (33) 



Mbar,thir) = 



4ai 



UrUpG 
ad$c , 



feT(r) ^ / dlnpgasjr) _^ d\nT{r) 



dlnr 



dlnr 



+ 



4ai 2a"±'c / s 
versus the observed mass contributions 

Mbar,obs{r) = Mgasir) + Mgal{r) + McDgal{r) 



(40) 



(41) 



Since not all the data involved in the above estimate have 
measurable errors, we cannot perform an exact x-square 
minimization: Actually, we can minimize the quantity: 



2 

X = 



{Mbar,obs — Mbar,theo)^ 



1 ^ 

i=l 



Mb ar,the 



(42) 



where N is the number of data and Hp = 2 the free pa- 
rameters of the model. We minimize the x-square using 
the Markov Chain Monte Carlo Method (MCMC). For each 
cluster, we have run various chains to set the best param- 
eters of the used algorithm, the Metropolis-Hastings one: 
starting from an initial parameter vector p (in our case 
p = (ai, 02)), we generate a new trial point p' from a tested 
proposal density g(p',p), which represents the conditional 
probability to get p', given p. This new point is accepted 
with probability 
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Table 2. Column 1; Cluster name. Column 2: first derivative coefficient, ai, of f(R) series. ColumnS: la confidence 
interval for ai. Column 4: second derivative coefficient, 02, of f(R) series. Column 5: la confidence interval for 02. 
Column 6: characteristic length, L, of the modified gravitational potential, derived from ai and 02- Column 7 : la 
confidence interval for L. 



name 


ai 


[ai - la, 


ai + la] 


02 


[02 — la, a2 + la] 


L 


[L -la, L + la] 










(kpc2) 


(kpc2) 


(kpc) 


(kpc) 


A133 


0.085 


[0.078, 


0.091] 


-4.98 • lO'"* 


[-2.38 • lO", -1.38 ■ 10^] 


591.78 


[323.34, 1259.50] 


A262 


0.065 


[0.061, 


0.071] 


-10.63 


[-57.65, -3.17] 


31.40 


[17.28, 71.10] 


A383 


0.099 


[0.093, 


0.108] 


-9.01 • 10^ 


[-4.10 • 10^, -3.14 • 102] 


234.13 


[142.10, 478.06] 


A478 


0.117 


[0.114, 


0.122] 


-4.61 • 10^ 


[-1.01 • 10"*, -2.51 • 10^] 


484.83 


[363.29, 707.73] 


A907 


0.129 


[0.125, 


0.136] 


-5.77 • 10^ 


[-1.54-104, -2.83-103] 


517.30 


[368.84, 825.00] 


A1413 


0.115 


[0.110, 


0.119] 


-9.45 • 10* 


[-4.26- IQB, -3.46- 10*] 


2224.57 


[1365.40, 4681.21] 


A1795 


0.093 


[0.084, 


0.103] 


-1.54 • 10^ 


[-1.01 • 10*, -2.49 • 10^] 


315.44 


[133.31, 769.17] 


A1991 


0.074 


[0.072, 


0.081] 


-50.69 


[-3.42 ■ 102, -1.3] 


64.00 


[32.63, 159.40] 


A2029 


0.129 


[0.123, 


0.134] 


-2.10 ■ 10"' 


[-7.95 ■ 10*, -8.44 ■ 10'^] 


988.85 


[637.71, 1890.07] 


A2390 


0.149 


[0.146, 


0.152] 


-1.40 ■ 108 


[-5.71 ■ 10**, -4.46 ■ 10^] 


7490.80 


[4245.74, 15715.60] 


MKW4 


0.054 


[0.049, 


0.060] 


-23.63 


[-1.15 • 10^, -8.13] 


51.31 


[30.44, 110.68] 


RXJ1159 


0.048 


[0.047, 


0.052] 


-18.33 


[-1.35 - 10^, -4.18] 


47.72 


[22.86, 125.96] 



f . /, L(d|p')P(p')g(p',p) 

1 L{d\p)P{p)q{p,p') 



where d are the data, L(d|p') oc exp(— x^/2) is the likeh- 
hood function, P(p) is the prior on the parameters. In our 
case, the prior on the fit parameters is related to Eq. (13): 
being L a length, we need to force the ratio axja^ to be 
positive. The proposal density is Gaussian symmetric with 
respect of the two vectors p and p', namely 9(p, p') oc 
exp(— Ap^/2(T^), with Ap = p — p'; we decide to fix the 
dispersion a of any trial distribution of parameters equal 
to 20% of trial a\ and ai at any step. This means that 
the parameter a reduces to the ratio between the likelihood 
functions. We have run one chain of 10® points for every 
cluster; the convergence of the chains has been tested using 
the power spectrum analysis from (Dunklcy 2005). The key 
idea of this method is, at the same time, simple and pow- 
erful: if we take the "power spectra of the MCMC samples, 
we will have a great correlation on small scales but, when 
the chain reaches convergence, the spectrum becomes fiat 
(like a white noise spectrum); so that, by checking the spec- 
trum of just one chain (instead of many parallel chains as in 
Gelmann-Rubin test) will be sufficient to assess the reached 
convergence. Remanding to (Dunklcy 2005) for a detailed 
discussion of all the mathematical steps. Here we calculate 
the discrete power spectrum of the chains: 



j |2 



with 



1 ^""^ 
^ n 



exp 



.277 7 



(43) 



(44) 



where A'^ and Xn are the length and the element of the sam- 
ple from the MCMC, respectively, j = 1, . . . , ^ - 1. The 
wavenumber kj of the spectrum is related to the index j 
by the relation kj = 2^ . Then we fit it with the analytical 
template: 



P(t)-p. 

or in the equivalent logarithmic form: 



(45) 



In Pj = In Po + In 



{k*/k, 



1 + ik*/k,)^ 



-l + rj 



(46) 



where 7 = 0.57216 is the Euler-Mascheroni number and 
rj are random measurement errors with < rj and 
< riVj >= Sij-K^ /6. From the fit, we estimate the two fun- 
damental parameters, Po and j* (the index corresponding to 
fc*). The first one is the value of the power spectrum extrap- 
olated for fe — > and, from it, we can derive the convergence 

ratio from r « --^; if r < 0.01, we can assume that the con- 
TV 

vcrgcncc is reached. The second parameter is related to the 
turning point from a power-law to a flat spectrum. It has 
to be > 20 in order to be sure that the number of points in 
the sample, coming from the convergence region, are more 
than the noise points. If these two conditions are verified for 
all the parameters, then the chain has reached the conver- 
gence and the statistics derived from MCMC well describes 
the underlying probability distribution (typical results are 
shown in Figs. 3, 4, 5. Following (Dunklcy 2005) prescrip- 
tions, we perform the fit over the range 1 < j < jmax , with 
jmax ~ Wj* , where a first estimation of j* can be obtained 
from a fit with jmax = 1000, and then performing a second 
iteration in order to have a better estimation of it. Even if 
the convergence is achieved after few thousand steps of the 
chain, we have decided to run longer chains of 10® points 
to reduce the noise from the histograms and avoid under- 
or over- estimations of errors on the parameters. The i — a 
confidence levels are easily estimated deriving them from the 
final sample the 15.87-th and 84.13-th quantiles (which de- 
fine the 68% confidence interval) for i = 1, the 2.28-th and 
97. 72-th quantiles (which define the 95% confidence inter- 
val) for i — 2 and the 0.13-th and 99.87-th quantiles (which 
define the 99% confidence interval) for i = 3. 

After the description of the method, let us now com- 
ment on the achieved results. 



6 RESULTS 

The numerical results of our fitting analysis are summa- 
rized in Table 2; we give the best fit values of the indepen- 
dent fitting parameters oi and 02, and of the gravitational 
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length L, considered as a function of the previous two quan- 
tities. In Figs. 3- 5, we give the typical results of fitting, with 
histograms and power spectrum of samples derived by the 
MCMC, to assess the reached convergence (fiat spectrum at 
large scales). 

The goodness and the properties of the fits are shown in 
Figs. 6- 17. The main property of our results is the presence 
of a typical scale for each cluster above which our model 
works really well (typical relative differences are less than 
5%), while for lower scale there is a groat difference. It is 
possible to see, by a rapid inspection, that this turning-point 
is located at a radius w 150 kpc. Except for extremely rich 
clusters, it is clear that this value is independent of the clus- 
ter, being similar for all clusters in our sample. 

There are two main independent explanations that 
could justify this trend: limits due to a break in the state of 
hydrostatic equilibrium or limits in the series expansion of 
the /(7?)-models. 

If the hypothesis of hydrostatic equilibrium is not cor- 
rect, then we are in a regime whore the fundamental relations 
Eqs. (23)- (28), are not working. As discussed in (Vikhlinin 
et al.2005), the central (70 kpc) region of most clusters is 
strongly affected by radiative cooling and thus its physical 
properties caimot directly be related to the depth of the 
cluster potential well. This means that, in this region, the 
gas is not in hydrostatic equilibrium but in a multi-phase 
state. In this case, the gas temperature cannot be used as a 
good standard tracer. 

We have also to consider another limit of our modelling: 
the requirement that the /(i?)-function is Taylor expand- 
able. The corrected gravitational potential which we have 
considered is derived in the weak field limit, which means 



R - Rq « 



ai 



(47) 



where Ro is the background value of the curvature. If this 
condition is not satisfied, the approach does not work (see 
(Capozziello, Stabile and Troisi 2008a) for a detailed discus- 
sion of this point). Considering that ai/a2 has the dimension 
of length'^ this condition defines the length scale where our 
series approximation can work. In other words, this indicates 
the limit in which the model can be compared with data. 

For the considered sample, the fit of the parameters ai 
and 02, spans the length range {19; 200} kpc (except for the 
richest clusters). It is evident that every galaxy cluster has 
a proper gravitational length scale. It is worth noticing that 
a similar situation, but at completely different scales, has 
been found out for low surface brightness galaxies modelled 
by /(i?)-gravity (Capozziello, Cardone and Troisi 2007). 

Considering the data at our disposal and the analy- 
sis which we have performed, it is not possible to quantify 
exactly the quantitative amount of these two different phe- 
nomena (i.e. the radiative cooling and the validity of the 
weak field limit). However, they are not mutually exclusive 
but should be considered in details in view of a more refined 
modelling ^. 



t Other secondary phenomena as cooling flows, merger and asym- 
metric shapes have to be considered in view of a detailed mod- 
elling of clusters. However, in this work, we are only interested to 
show that extended gravity could be a valid alternative to dark 
matter in order to explain the cluster dynamics. 




Pgas ( kg/m' ) 

Figure 2. Density vs ai: predictions on the behavior of ai. 

The horizontal black bold Une indicates the Newtonian-limit, 
ai — > 3/4 which we expect to be realized on scales comparable 
with Solar System. Vertical lines indicate typical approximated 
values of matter density (without dark matter) for different grav- 
itational structures; universe (large dashed) with critical den- 
sity Pcrit ^ 10~^® kg/m^; galaxy clusters (short dashed) with 
Pci ^ 10~^^ kg/m^; galaxies (dot-dashed) with pgg^i fs 10~^^ 



kg/m^; sun (dotted) with ps 
show the predicted trend for ai. 



10 kg/m . Arrows and boxes 



Similar issues are present also in (Brownstein 2006B): 
they use the the Metric - Skew - Tensor - Gravity (MSTG) 
as a generalization of the Einstein General Relativity and 
derive the gas mass profile of a sample of clusters with gas 
being the only baryonic component of the clusters. They 
consider some clusters included in our sample (in particu- 
lar, A133, A262, A478, A1413, A1795, A2029, MKW4) and 
they find the same different trend for r < 200 kpc, even if 
with a different behavior with respect to us: our model gives 
lower values than X-ray gas mass data while their model 
gives higher values with respect to X-ray gas mass data. 
This stresses the need for a more accurate modelling of the 
gravitational potential. 

However, our goal is to show that potential (12) is suit- 
able to fit the mass profile of galaxy clusters and that it 
comes from a self-consistent theory. 

In general, it can be shown that the weak field limit 
of extended theories of gravity has Yukawa-like correc- 
tions (Stelle 1978; Kenmoku, Okamoto and Shigemotol993). 
Specifically, given theory of gravity of order (2n + 2), 
the Yukawa corrections to the Newtonian potential are n 
(Quandt 1991). This means that if the effective Lagrangian 
of the theory is 



C = /(i?, aR,..u''R, 
we have 



"R)y 



'^(^) = — 



1 + y^, Qfcg 



(48) 



(49) 



Standard General Relativity, where Yukawa corrections are 
not present, is recovered for n = (second order theory) 
while the /(i?)-gravity is obtained for n = 1 (fourth-order 
theory). Any □ operator introduces two further derivation 
orders in the field equations. This kind of Lagrangian comes 
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out when quantum field theory is formulated on curved 
spacctime (Birrell & Daviesl982). In the series (49), G is 
the value of the gravitational constant considered at infin- 
ity, Lk is the interaction length of the fe-th component of 
the non-Newtonian corrections. The amplitude at of each 
component is normalized to the standard Newtonian term; 
the sign of at tells us if the corrections are attractive or re- 
pulsive (see (Willl993) for details). Moreover, the variation 
of the gravitational coupling is involved. In our case, we are 
taking into account only the first term of the series. It is the 
the leading term. Let us rewrite (12) as 



0(r) 



GM r 

1 + aie 

r 



r/Li 



(50) 



The effect of non-Newtonian term can be parameterized by 
{ai, Li} which could be a useful paramctcrisation which 
respect to our previous {ai, 02} or {Geff, L} with Geff = 
3G/(4oi). For large distances, where r ^ Li, the exponen- 
tial term vanishes and the gravitational coupling is G. If 
r <^ Li, the exponential becomes 1 and, by differentiating 
Eq.(50) and comparing with the gravitational force mea- 
sured in laboratory, we get 



Glab = G 



l + ai 



(-i 



-r/Li 



:G(l-hQ;i), 



(51) 



where Giab = 6.67 x 10 g cm s is the usual New- 
ton constant measured by Cavendish-like experiments. Of 
course, G and Giab coincide in the standaxd Newtonian grav- 
ity. It is worth noticing that, asymptotically, the inverse 
square law holds but the measured coupling constant differs 
by a factor (1 -|- ai). In general, any correction introduces a 
characteristic length that acts at a certain scale for the self- 
gravitating systems as in the case of galaxy cluster which we 
are examining here. The range of Lk of the fcth-component 
of non-Newtonian force can be identified with the mass nik 
of a pseudo-particle whose effective Compton's length can 
be defined as 

Lk = —. (52) 

mfeC 

The interpretation of this fact is that, in the weak energy 
limit, fundamental theories which attempt to unify gravity 
with the other forces introduce, in addition to the massless 
gravitori, particles with mass which also carry the gravita- 
tional interaction (Gibbons & Whitingl981). See, in particu- 
lar, (Capozziello et al. 2008b) for /(7?)-gravity. These masses 
are related to effective length scales which can be parame- 
terized as 

Lfc = 2 X 10-"^ f— cm. (53) 
\ ruk J 

There have been several attempts to experimentally con- 
strain Lk and Ok (and then irik) by experiments on scales 
in the range 1 cm < r < 1000 km, using different tech- 
niques (Fischbach et al. 1986; Speake & Quinnl988; Eck- 
hardtl993). In this case, the expected masses of particles 
which should carry the additional gravitational force are in 
the range 10~'^^eV < mk < 10~^eV. The general outcome 
of these experiments, even retaining only the term fc = 1, 
is that geophysical window between the laboratory and the 
astronomical scales has to be taken into account. In fact, the 
range 



is not excluded at all in this window. An interesting sug- 
gestion has been given by Fujii (Fujiil988), which proposed 
that the exponential deviation from the Newtonian standard 
potential could arise from the microscopic interaction which 
couples the nuclear isospin and the baryon number. 

The astrophysical counterparts of these non-Newtonian 
corrections seemed ruled out till some years ago due to the 
fact that experimental tests of General Relativity seemed to 
predict the Newtonian potential in the weak energy limit, 
"inside" the Solar System. However, as it has been shown, 
several alternative theories seem to evade the Solar System 
constraints (see (Capozziello et al.2008b) and the reference 
therein for recent results) and, furthermore, indications of an 
anomalous, long range acceleration revealed from the data 
analysis of Pioneer 10/11, Galileo, and Ulysses spacecrafts 
(which axe now almost outside the Solar System) mates 
these Yukawa-like corrections come again into play (An- 
derson et al.l998). Besides, it is possible to reproduce phc- 
nomenologically the flat rotation curves of spiral galaxies 
considering the values 



ai = -0.92 , Li ~ 40 kpc . 



(55) 



10" 



Li 



l(f ^ 10^ m , 



(54) 



The main hypothesis of this approach is that the additional 
gravitational interaction is carried by some ultra-soft boson 
whose range of mass is mi ~ lO"^'^ -i- 10~^*eV. The action 
of this boson becomes efficient at galactic scales without the 
request of enormous amounts of dark matter to stabilize the 
systems (Sandersl990). 

Furthermore, it is possible to use a combination of two 
exponential correction terms and give a detailed explanar 
tion of the kinematics of galaxies and galaxy clusters, again 
without dark matter model (Eckliardtl993). 

It is worthwhile to note that both the spacecrafts mea- 
surements and galactic rotation curves indications come 
from "outside" the usual Solar System boundaries used up 
to now to test General Relativity. However, the above re- 
sults do not come from any fundamental theory to explain 
the outcome of Yukawa corrections. In their contexts, these 
terms axe phenomenological. 

Another important remark in this direction deserves the 
fact that some authors (McGaugh2000) interpret also the ex- 
periments on cosmic microwave background like the experi- 
ment BOOMERANG and WMAP (de Bernaxdis et al. 2000; 
Spergel et al. 2003) in the framework of modified Newtonian 
dynamics again without invoking any dark matter model. 

All these facts point towards the line of thinking that 
also corrections to the standard gravity have to be seriously 
taken into account beside dark matter searches. 

In our case, the parameters ai_2, which determine 
the gravitational correction and the gravitational coupling, 
come out "directly" from a field theory with the only re- 
quirement that the effective action of gravity could be more 
general than the Hilbert-Einstein theory f{R) = R. This 
main hypothesis comes from fundamental physics motiva- 
tions due to the fact that any unification scheme or quan- 
tum field theory on curved space have to take into ac- 
count higher order terms in curvature invariants (Birrell & 
Daviesl982). Besides, several recent results point out that 
such corrections have a main role also at astrophysical and 
cosmological scales. For a detailed discussion, see (Nojiri & 
Odintsov 2007; Capozziello & Francaviglia 2008; Sotiriou & 
Faraoni2008). 
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With this philosophy in mind, we have plotted the trend 
of ai as a function of the density in Fig.2. As one can see, 
its values are strongly constrained in a narrow region of the 
parameter space, so that oi can be considered a "tracer" for 
the size of gravitational structures. The value of ai range 
between {0.8 0.12} for larger clusters and {0.4 0.6} 
for poorer structures (i.e. galaxy groups like MKW4 and 
RXJ1159). We expect a particular trend when applying the 
model to different gravitational structures. In Fig. 2, wc give 
characteristic values of density which range from the biggest 
structure, the observed Universe (large dashed vertical line), 
to the smallest one, the Sun (vertical dotted line), through 
intermediate steps like clusters (vertical short dashed line) 
and galaxies (vertical dot-dashed line). The bold black hori- 
zontal line represents the Newtonian limit ai =3/4 and the 
boxes indicate the possible values of ai that we obtain by 
applying our theoretical model to different structures. 

Similar considerations hold also for the characteristic 
gravitational length L directly related to both ai and 02. 
The parameter 02 shows a very large range of variation 
{—10^ -. — 10} with respect to the density (and the mass) of 
the clusters. The value of L changes with the sizes of grav- 
itational structure (see Fig. 18), so it can be considered, 
beside the Schwarzschild radius, a sort of additional grav- 
itational radius. Particular care must be taken when con- 
sidering Abell 2390, which shows large cavities in the X-ray 
surface brightness distribution, and whose central region, 
highly asymmetric, is not expected to be in hydrostatic equi- 
librium. All results at small and medium radii for this cluster 
could hence be strongly biased by these effects (Vikhlinin et 
al.2006); the same will hold for the resulting exceptionally 
high value of L. Fig. 18 shows how observational proper- 
ties of the cluster, which well characterize its gravitational 
potential (such as the average temperature and the total 
cluster mass within rsoo, plotted in the left and right panel, 
respectively), well correlate with the characteristic gravitar- 
tional length L. 

For clusters, we can define a gas-density- weighted and a 
gas-mass-weighted mean, both depending on the series pa- 
rameters ai,2- We have: 

< L>p = 318 kpc < a2 >p= -3.40 • 10* 
<L>M = 2738 kpc < a2 >M= -4.15 • 10^ (56) 

It is straightforward to note the correlation with the sizes 
of the cluster cD-dominated-central region and the "grav- 
itational" interaction length of the whole cluster. In other 
words, the parameters ai,2, directly related to the first and 
second derivative of a given analytic /(i?)-model determine 
the characteristic sizes of the self gravitating structures. 



7 DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the possibility that the 
high observational mass-to-light ratio of galaxy clusters 
could be addressed by f{R)- gravity without assuming huge 
amounts of dark matter. We point out that this proposal 
comes out from the fact that, up to now, no definitive can- 
didate for dark-matter has been observed at fundamental 
level and then alternative solutions to the problem should 
be viable. Furthermore, several results in /(i?)-gravity seem 
to confirm that valid alternatives to ACDM can be achieved 



in cosmology. Besides, as discussed in the Introduction, the 
rotation curves of spiral galaxies can be explained in the 
weak field limit of /(i?)-gravity. Results of our analysis go 
in this direction. 

We have chosen a sample of relaxed galaxy clusters 
for which accurate spectroscopic temperature measurements 
and gas mass profiles are available. For the sake of simplic- 
ity, and considered the sample at our disposal, every clus- 
ter has been modelled as a self bound gravitational sys- 
tem with spherical symmetry and in hydrostatic equilib- 
rium. The mass distribution has been described by a cor- 
rected gravitational potential obtained from a generic ana/- 
lytic /(ii)-theory. In fact, as soon as f{R) ^ R, Yukawa- 
like exponential corrections emerge in the weak field limit 
while the standard Newtonian potential is recovered only 
for f{R) = R, the Hilbert-Einstein theory. 

Our goal has been to analyze if the dark-matter content 

of clusters can be addressed by these correction potential 
terms. As discussed in detail in the previous section and how 
it is possible to see by a rapid inspection of Figs. 6- 17, all 
the clusters of the sample axe consistent with the proposed 
model at Icr confidence level. This shows, at least qualita- 
tively, that the high mass-to-light ratio of clusters can be 
explained by using a modified gravitational potential. The 
good agreement is achieved on distance scales starting from 
150 kpc up to 1000 kpc. The differences observed at smaller 
scales can be ascribed to non-gravitational phenomena, such 
as cooling flows, or to the fact that the gas mass is not a 
good tracer at this scales. The remarkable result is that we 
have obtained a consistent agreement with data only using 
the corrected gravitational potential in a large range of radii. 
In order to put in evidence this trend, we have plotted the 
baryonic mass vs radii considering, for each cluster, the scale 
where the trend is clearly evident. 

In our knowledge, the fact that /(il)-gravity could work 
at these scales has been only supposed but never achieved 
by a direct fltting with data (see (Lobo 2008) for a review) . 
Starting from the series coefficients ai and a2! it is possi- 
ble to state that, at cluster scales, two characteristic sizes 
emerge from the weak field limit of the theory. However, at 
smaller scales, e.g. Solar System scales, standard Newtonian 
gravity has to be dominant in agreement with observations. 

In conclusion, if our considerations are right, gravita- 
tional interaction depends on the scale and the infrared limit 
is led by the series coefficient of the considered effective 
gravitational Lagrangian. Roughly speaking, we expect that 
starting from cluster scale to galaxy scale, and then down 
to smaller scales as Solar System or Earth, the terms of the 
series lead the clustering of self-gravitating systems beside 
other non-gravitational phenomena. In our case, the New- 
tonian limit is recovered for ai 3/4 and L(ai,a2) ^ r 
at small scales and for L{ai,a2) <C r at large scales. In the 
first case, the gravitational coupling has to be redifined, in 
the second Goo — G- In these limits, the linear Ricci term 
is dominant in the gravitational Lagrangian and the New- 
tonian gravity is restored (Quandt 1991). Reversing the ar- 
gument, this could be the starting point to achieve a theory 
capable of explaining the strong segregation in masses and 
sizes of gravitationally-bound systems. 
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Figure 3. Left panel: histogram of the sample points for parameter ai in Abell 383 coming out the MCMC implementation used to 
estimate best fit values and errors for our fitting procedure as described in § 5.5. Binning (horizontal axis) and relative frequencies 
(vertical axis) are given by automatic procedure from Mathematica6.0. Right panel: power spectrum test on sample chain for parameter 
ai using the method described in § 5.5. Black line is the logarithm of the analytical template Eq. (45) for power spectrum; gray line is 
the discrete power spectrum obtained using Eq. (43) - (44). 
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Figure 4. Abell 383: histogram (left) and power spectrum test (right) on sample chain for parameter 02. 




Figure 5. Abell 383: histogram (left) and power spectrum test (right) on sample chain for parameter L. 
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Figure 6. Baryonic mass vs radii for Abell A133. Dashed line is the experimental-observed estimation Eq. (41) of baryonic matter 
component (i.e. gas, galaxies and cD-galaxy); solid line is the theoretical estimation Eq. (40) for baryonic matter component. Dotted 
lines are the 1-<t confidence levels given by errors on fitting parameters in the left panel; and from fitting parameter plus statistical errors 
on mass profiles as discussed in § 5.4 in the right panel. 
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Figure 7. Same of Fig. 6 but for cluster Abell 262. 
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Figure 8. Same of Fig. 6 but for cluster Abell 383. 
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Figure 9. Same of Fig. 6 but for cluster Abell 478. 
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Figure 10. Same of Fig. 6 but for cluster Abell 907. 
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Figure 11. Same of Fig. 6 but for cluster Abell 1413. 
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Figure 12. Same of Fig. 6 but for cluster Abell 1795. 
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Figure 14. Same of Fig. 6 but for cluster Abell 2029. 
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Figure 15. Same of Fig. 6 but for cluster Abell 2390. 
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Figure 18. Single temperature fit to the total cluster spectrum (left panel) and total cluster mass within rsoo (given as a function of 
Mq) (right panel) are plotted as a function of the characteristic gravitational length L. Temperature and mass values are from Vikhlinin 
et al.2006. 
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